Infinite families of (quasi)-Hermitian Hamiltonians associated with 
exceptional X m Jacobi polynomials 



a* 

00 



Bikashkali Midya* and Barnana Royt 

Physics and Applied Mathematics Unit, 
Indian Statistical Institute, 
Kolkata 700108, India. 



Abstract 

Using an appropriate change of variable, the Schrodinger equation is transformed into a second-order 



r~| [ differential equation satisfied by recently discovered Jacobi type X m exceptional orthogonal polynomials. 

This facilitates the derivation of infinite families of exactly solvable Hermitian as well as non-Hermitian 
trigonometric and hyperbolic Scarf potentials whose bound state solutions are associated with the aforesaid 
exceptional orthogonal polynomials. These infinite families of potentials are shown to be extensions of 
the conventional trigonometric and hyperbolic Scarf potentials by the addition of some rational terms 
characterized by the presence of classical Jacobi polynomials. All the members of a particular family of 
these 'rationally extended polynomial-dependent' potentials have the same energy spectrum and possess 
translational shape invariant symmetry. The obtained non-Hermitian potentials are shown to be quasi- 
Hermitian in nature ensuring the reality of the associated energy spectra. 
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I. INTRODUCTION 



In recent years there is a surge of interest in the study of exactly solvable quantum systems 
associated with exceptional orthogonal polynomials introduced in the seminal paper by Gomez- 
Ullate et al [1, 2]. The exceptional X m orthogonal polynomials (EOP) are the solutions of 
second-order Sturm-Liouville eigenvalue problem with rational coefficients. A distinguishing 
property of these polynomials is that the lowest eigen polynomial of the sequence need not be 
of degree zero, even though the full set of eigenfunctions still forms a basis of the weighted C? space. 

Exceptional X\ Laguerre or Jacobi type polynomials were shown to be eigenfunctions of ratio- 
nally extended radial oscillator or Scarf I potentials by the point canonical transformation method 
[3] and by supersymmetric quantum mechanical method (SUSYQM) [4, 5]. The latter is based on 
a reparametrization of the starting Hamiltonian and redefining the conventional superpotential in 
terms of modified couplings, while allowing for the presence of extra rational contribution expressed 
in terms of some polynomial function. Depending on the degree of the latter, the wave functions 
of the extended potentials are shown to be X\ or X 2 Laguerre, Jacobi exceptional orthogonal 
polynomials. Employing similar approach in the context of second order supersymmetry, rational 
extensions of the generalized Poschl-Teller, Scarf I potentials with the corresponding bound state 
wavefunctions given in terms of exceptional orthogonal polynomials have been constructed [5]. 
Construction of two distinct families of Laguerre and Jacobi type X m , m = 1,2,3..., exceptional 
orthogonal polynomials as eigenfunctions of infinitely many shape invariant potentials by deform- 
ing the radial oscillator and the Hyperbolic/Trigonometric Poschl Teller potentials was made in 
refs.[6, 7] and the properties of these X m exceptional orthogonal polynomials have been thoroughly 
studied subsequently [8-12]. The multi-indexed version X mi m2j mfe of exceptional X m orthog- 
onal polynomials are constructed by using multi-step Darboux transformation [13], Crum-Adler 
mechanism [14], higher-order SUSYQM [15] and multi-step Darboux-Backlund transformation [16]. 

X m exceptional orthogonal polynomials were shown to be obtainable through several other 
approaches such as Darboux-Crum transformation [17-19], Darboux-Backlund transformation 
[16, 21, 22] and prepotential approach [23, 24]. In fact, X m Laguerre polynomials have been 
characterized in terms of an isospectral Darboux transformation. The shape invariance of these 
polynomial families is shown to be a direct consequence of the permutability property of the 
Darboux transformation [17]. It is shown to be possible to obtain rational extensions of every 
translationally shape invariant potentials via Darboux Backlund transformations based on negative 
eigenfunctions built from excited states of the initial Hamiltonian. The regularity of corresponding 
Riccati-Schrodinger functions can be verified by combining disconjugacy theorems and asymptotic 
analysis of the wave functions. The potentials obtained by this method are either strictly 
isospectral or quasi isospectral extensions of the initial ones and have enlarged shape invariance 
property. All the quantal systems related to the exceptional Laguerre and Jacobi polynomials can 



2 



be constructed by the prepotential approach [23, 24] which does not need the concept of shape 
invariance and the Darboux-Crum transformation. The prepotential, the deforming function, the 
potential, the eigenfunctions and eigenvalues can be derived within the same framework. The 
exceptional polynomials obtained in this way are expressible as a bilinear combination of the 
deforming functions and its derivatives. 

Recently, EOP have been studied in diverse scenarios. For instance, these new polynomials 
were shown to be associated with the solutions of some conditionally exactly solvable potentials 
[25, 26], solutions for position dependent mass systems [27], main part of the eigenfunctions 
of the Dirac equations coupled minimally and non-minimally with some external fields and the 
Fokker-Planck equations [28]. Four types of infinitely many exactly solvable Fokker-Planck 
equation which are the generalized, or deformed versions of the Rayleigh process and the Jacobi 
process are also shown to related to EOP [29]. The Structure of the X m Laguerre polynomials 
was considered within the quantum Hamilton-Jacobi formalism [30], within N-fold supersymmetry 
[31] and its dynamical breaking in the context of position dependent mass scenario [32]. 

On the other hand, relaxation of Hermiticity for the reality of the discrete spectrum of a 
quantum mechanical Hamiltonian has given rise to some very interesting investigations in the last 
few years stepped up by a conjecture of Bender and Boettcher [33] that PT-symmetric non- 
Hermitian Hamiltonians could possess real bound state spectrum. The complex PT-sym metric 
Hamiltonians possess a real discrete spectrum if the energy eigenstates are also the eigenstates of 
VT; otherwise, the VT symmetry is said to be spontaneously broken and the energy eigenvalues 
occur in complex conjugate pairs [34]. Here V stands for parity transformation and T stands 
for time reversal. Subsequently, Mostafazadeh [35] showed that the reality of the spectrum of 
a VT symmetric Hamiltonian is ensured if the Hamiltonian H is Hermitian with respect to a 
positive definite inner product (.\r)\.) on the Hilbert space "H in which H is acting. This renders 
the Hamiltonian H to be pseudo-Hermitian [36] = r]Hr]~ l , where the Hermitian linear 
automorphism r] : H — > 7i is bounded and positive definite. Another equivalent condition for 
the reality of the energy spectrum of H is the quasi-Hermiticity [37-39], i.e. the existence of 
a invertible operator p such that H = pHp^ 1 is Hermitian with respect to usual inner product 
(.,.). Quasi-Hermitian Hamiltonian shares the same energy spectrum of the equivalent Hermitian 
Hamiltonian H and the wave functions are obtained by operating on those of H . Most of the 
analytically solvable non-Hermitian Hamiltonians are constructed by either making the coupling 
constant of the potential present in the relevant Schrodinger Hamiltonian imaginary [40, 41] 
or by shifting the coordinate with an imaginary constant [42, 43]. Several of these classes of 
Hamiltonians are argued to be pseudo-Hermitian under 77 = e~ ap [44]. For a real a and p = —i-^, 
the operator rj shifts the coordinate x to x + ia. 
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In this paper we have found infinite families of exactly solvable Hermitian as well as non- 
Hermitian Hamiltonians whose bound state wave functions are given in terms of exceptional X m 
Jacobi orthogonal polynomials. The non-Hermitian Hamiltonians, obtained here, are shown to 
be quasi-Hermitian in nature rendering these Hamiltonians to have real eigenvalues. Also, using 
SUSYQM, we have shown that all the infinite family of Hamiltonians have translational shape 
invariant property. The motivation for doing this comes from the fact that compared to the 
exactly solvable real potentials associated with exceptional orthogonal X m polynomials, a little 
attention [4, 45] has been given to look for complex potentials having these polynomials as their 
bound state eigenfunctions. In fact a few complex potentials have been found involving only X\ 
type Jacobi, Laguerre exceptional orthogonal polynomials. Considering the recent flurry of interest 
in non-Hermitian Hamiltonians with VT symmetry and/or pseudo-Hermiticity the present study 
deserves to be worthwhile. The organization of the paper is as follows: In section II, some properties 
of X m Jacobi EOPs, which are necessary to analyze the results, are given. Construction of infinite 
families of Hermitian and non-Hermitian Hamiltonians by transforming the relevant Schrodinger 
equation into the differential equation satisfied by the X m Jacobi polynomials are done in section 
III. In section IV, the shape invariance properties of the obtained potentials are discussed within the 
framework of supersymmetric quantum mechanics. Finally, section V is devoted to the summary 
and outlook. 



II. EXCEPTIONAL X m JACOBI ORTHOGONAL POLYNOMIALS 

We present here some properties of X m Jacobi polynomials which are studied in detail in ref. [19]. 
For a fixed integer m > 1 and real a,b> — 1, the exceptional X m Jacobi orthogonal polynomials 
Pn ,b ' m \x), n — m,m + l,m + 2..., satisfies the differential equation 

{l - x2)d ^ + Ql{x) tx +Rl{x)v = ^ {21) 



where Qi(x) and R\(x) can be written in terms of classical Jacobi orthogonal polynomial Pm' b \x), 
as 

Q 1 (x) = (a-b-m + l) 1 a \l n -} ] -(fl + i)(i+ x ) + (6+i)(i- x ) 

Pm ' [X) 

(2.2) 

(1 - x)P ( -~ a,b ^(x) 

RJx) = b(a-b-m + l) - — , ™~\ + n 2 + n(a + b - 2m + 1) - 2bm. 

P^-^ix) 

For n,l > m, this new family of polynomials are orthonormal with respect to the weight function 
W a ,b,m(x) = ^Za-S-i^ ■ Evidently, for £ 2 -orthonormality to hold, the denominator of this weight 
function should be non-zero in — 1 < x < 1. According to the references [19, 46] the polynomial 
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p(o+l,6-l) 



x) has no zeroes in [—1, 1] if and only if the following conditions hold simultaneously: 

(2.3) 



(i) b ^ 0, a, a — b — m + 1 ^ {0, 1, m — 1}, 

(ii) a > m — 2, sgn(a — m + 1) = sgn(6), 



where sgn(x) is the signum function. With these restrictions the £ 2 -norms of the exceptional X„ 
Jacobi polynomials are given by 



1 (1 -x) a (l +x) 



P 



(-a-1,6-1) 



X) 



(2.4) 



2 a+b+1 (n + b){n - 2m + a + l)T(n - m + a + 2)T(n + b) 
(2n — 2m + a + b + l)(n — m + a + l) 2 (n — m)\ T{n — m + a + b+1) 



Jnh 



where 5 n i is the Kronecker delta function. Moreover, the X m Jacobi orthogonal polynomials 



y(a,b,m) 



x) are related to the classical Jacobi polynomials Pn ,b \x) by the following relations 



p(a,b,m) = (_-Qr 



{x _ l)p (-a-l, b -l)pia+ 



+2,< 



a+b+j+1 
2(a + j + l) - -'-™ 

■ a ~ m + 1 p(-a-2,b)p(a+l,&-l) 



(2.5a) 
(2.5b) 



j = n — m > 0. 



a+j+l~ m 

It is worth mentioning here that following identities of classical Jacobi polynomial are extensively 
used to simplify the results presented in the later part of the paper 

(1 - x 2 )(a + b + m + 2)P^' b+2) + 2[b-a-(a + b + 2)x]P^ h+1) + 4mP^ 6) = (2.6a) 
(x - l)(a + b + m + l)Pi a +i' 6+1) = 2(a + m)P^- 1 ' b+1 ^ - 2aP^ b) (2.6b) 
^ a ' 6 - 1} " Pt 1 * = P£l (2.6c) 



<7 r 



r(a + 6 + m + r + Z) (a+rih+r) 
2T(a + 6 + n + 1) ' 



(2.6d) 



The identity (2.6a), which is absent in the handbook [47], directly follows from the differential 
equation of the classical Jacobi polynomials Pm' b \x). 



III. INFINITE FAMILY OF HERMITIAN AND QUASI-HERMITIAN HAMILTONIANS ASSO- 
CIATED WITH JACOBI TYPE X m EOPS 

Here we shall derive some exactly solvable Hermitian as well as non-Hermitian potentials whose 
bound state solutions are connected to Jacobi type X m EOP. For this we have used an appropriate 
coordinate transformation which transform the Schrodinger equation into a differential equation of 
X m Jacobi EOP. In the following we first briefly describe this method. 
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A. Construction method 



Following the discussion of [45, 48], we consider transformation of the Schrodinger equation (with 

h = 2m = \) 

Hifj(x) = ~ d ^f + V(x)if>(x) = Ei/>{x). (3.1) 
into a second-order differential equation satisfied by a special function F(g). For this the solution 
ij)(x) of the equation (3.1) can be assumed as 

1>(x) ~ f(x) F{g(x)). (3.2) 

Substituting this form of ip(x) in (3.1), we obtain the following second-order ordinary differential 
equation of the special function F(g) 

d 2 F dF 

-^ + Q(g)^ + R(g)F(g) = 0, (3.3) 
where Q(g(x)) and R(g(x)) are given by 

m = (3-) 

respectively. After doing some algebraic manipulations, the equations (3.4a) and (3.4b) reduce to 

f{x)*g\x)-V*e yQ{a)d \ (3.5a) 
9'" 3 (g"\ 2 , /2 / ldQ 1 2 



E - v ^-l{7) +9 { R -2^-l Q )- (3 - 5b) 

The above expression of E — V(x) contains the unknown functions Q(g), R(g) and the change 
of variable g(x). A particular choice of the special function F(g) provides the complete functional 
forms of the first two unknowns, whereas a simple way of finding the unknown function g(x) 
has been first proposed by Bhattacharjee and Sudarshan [49]. According to them, if there is a 
constant (E) on the left-hand side of (3.5b), there must be one on the right-hand side too. Hence, 
a specific special function F(g) and a reasonable choice of g{x) make the Schrodinger equation 
(3.1) an exactly solvable problem for the potential V(x) (which is different for different choice of 
the special function). The solution of the corresponding Schrodinger equation can be obtained 
using the equations (3.2) and (3.5a) as 

^(x)^g'( x )-^e UQ(a)da F(g(x)). (3.6) 

This procedure has been exploited in ref.[48] and a systematic list of shape invariant real potentials 
were obtained whose solutions are related to the classical orthogonal polynomials. Analogous study, 
corresponding to X 1 EOPs, has been carried out in [3] for Hermitian systems and in [45] for quasi- 
Hermitian systems. In the following we use the above mentioned method to obtain infinitely many 
exactly solvable Hermitian as well as non-Hermitian Hamiltonians whose bound state wave functions 
are connected with exceptional X m Jacobi orthogonal polynomials for arbitrary m = 1,2,3.... 
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B. Infinite family of Hermitian Hamiltonians 



We choose the special function F{g) to be Jacobi type X m EOP e.g. F{g) ~ Pn (g). In this 
case we have from equations (2.1) and (3.3) 



Q( 9 ) = ^4, R( 9 )- Rl{9) 



(3.7) 



1-gV ^ 1-g^ 

where Qi(g) and R\(g) are given in equation (2.2). Using equations (2.2) and (3.7), equation 
(3.5b) reduces to 

1 - b 2 g' 2 



,2 g ,2 



~ V[X) - 27 " 4 [jj + 4 (l-g) 2+ 4 (1 + g) 2 



+ 



In 2 + 2n(a + b - 2m + 1) + 2m(a - 36 - m + 1) + (a + 1)(6 + 1) g'' 



1-9 2 



+ 



{a-b-m + l)[a + b+{a-b + l)g]g' 2 P^iHg) (a-b-m + l) 2 g 



Pt a i b) (9) 



1-9 2 Pk- a ~ hb ~ 1} (g) 2 

At this point we choose g' 2 /(l — g 2 ) = constant = ce I - {0}, which is satisfied by 

g(x) = sin(A/c x + d). 



p(-a-l,b-l) 



(3.8) 
(3.9) 



(9). 



Without loss of generality we first choose d = 0. Nonzero real values of d do not make any 
significant difference in the potential and its solutions. The nonzero complex values of d give rise 
to complex potentials with real spectra as will be shown later. The choice (3.9) of g(x) helps to 
generate a constant term in the right hand side corresponding to the energy in the left hand side 
of equation (3.8). Here two cases may arises, depending on the sign of the chosen constant c. In 
the following, we consider two cases c > and < separately, c > gives the infinite numbers 
of rationally extended trigonometric Scarf potential while c < corresponds the hyperbolic Scarf 
potential family. 



Case I : c > , extended trigonometric Scarf potential family 

Let c = k 2 , k 7^ 0. In this case, we have g(x) = sin kx. Substituting g(x) 
(3.8) and separating out the potential V(x) and energy E, we have 



sin/cx in equation 



V^ m \x) = k < " 2a +2b — — sec 2 kx — k ^ 2 a - sec kx tan kx — 1k 2 m{a — b — m + 1) 



k 2 (a — b — m + l)[a + b + (a — 6+1) sin kx] 



P, 



(-o,6), 



m—1 



sin kx) 



(-a-1,6-1) 



(sin kx) 



+ 



k 2 (a — b — m + l) 2 cos 2 kx 



P, 



(-0,6), 



m— 1 



sin kx) 



P 



(-a-1,6-1) 



sin kx) 



71 71 

2k <X< 2k 



(3.10a) 
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k 2 

#M = — (2n-2m + a + b + lf, n = m, m + 1, m + 2, ... (3.10b) 



Also, the expression of the wave function ip(x) given in equation (3.6) reduces to 

(1 — sin /cx) 5+3(1 + sin/cx)5 + 4 
Pm~ 1 ' b ~ 1 (sin fcx) 

The normalization constant can be obtained using equation (2.4), in the following form 



,-'(*) = * ~ Sm ™!X^;. ^ RX >' ' P^"'(3in fa), „> m (3.11) 



jy( m ) (2n — 2m + a + 6 + l)(n — m + a + l) 2 (n — m)l T(n — m + a + 6+1) 

n ' [ 2 a + b + 1 (n + b)(n-2m + a + l)T(n - m + a + 2)r(n + 6) 

(3.12) 

At this point few remarks on the obtained potentials V (jn \x) are worth to be mentioned. 

• The new potentials V^ m '{x), m = 0,1,2... become singular at the zeroes of the Jacobi 
polynomial Pm~ 1,b ~ (x) in the interval (— ;fr,5r)- These singularities can be avoided by 
restricting the potential parameters a,b according to the conditions mentioned in equation 
(2.3). 

• The potentials V (jn \x) given in (3.10a) are infinite in number because each integer values 
of m > gives rise to a new exactly solvable potential. 

• For m = the last three terms in the expression of V^ m \x) vanish and one is left with the well 
known trigonometric Scarf potential V^i^x) = k ^ 2a +2h sec 2 kx— k - ~ a ^ sec fcxtan fcx. 
It's solutions are given by classical Jacobi polynomials Pm' b \x) as can be easily verified using 
equations (3.11) and (2.5a). All other choices of m make the potentials V^ix) as the 
extensions of V <y0 \x) by addition of some rational terms which are explicitly depends on the 
classical Jacobi polynomials Pm' b \x). Therefore we can interpret V^ m \x) as the rationally 
extended polynomial-depended trigonometric Scarf potential family. 

• All the members V^(x), V^(x), V^(x)... of the polynomial-depended trigonometric Scarf 
potential family V <jn \x) are isospectral i.e. they share the same energy spectrum E^j = 
^-(2n' + a + b + l) 2 , n' = 0, 1, 2..., where we have reset the quantum number n = n' + m. 

In the following we are giving some simple cases of the results obtained in equations (3.10a), 
(3.10b) and (3.11). 

(a) m = : For m = 0, we have from equations (3.10a) and (3.10b) 



T/(0) . . k 2 (2a 2 + 2b 2 -l) 2l k 2 (b 2 -a 2 ) 

V K ' [x) = sec kx sec kx tan kx 



k 2 

E = —(2n + a + b+l) 2 
4 



(3.13) 
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By redefining the parameters a = a — (3 — 1/2, b = a+/3 — 1/2, the above equation (3.13) reduces 
to the well known trigonometric Scarf potential [48, 50] and its bound state energy spectrum given 
by 

V {0 \x) = k 2 [a(a - 1) + (3 2 ] sec 2 kx - k 2 (3(2a - 1) sec kxtan kx 

(3.14) 

= k 2 (n + a) 2 , n = 0, Ir- 
respectively. Corresponding wave functions can be written in terms of classical Jacobi polynomials 

= Af^(l-smkxY^\l + smkx)^ pi a '^' a+fi ^\sm kx) (3.15) 

(b) m = 1 : For m — 1 and the same choices of parameters mentioned in the case of m — 0, 
equation (3.8) reduces to 

V^(x) = k 2 [a(a - 1) + /3 2 ] sec 2 kx - k 2 /3(2a - 1) seckxt^kx + J^f^ - gfegggi 
= fc 2 (n + a - I) 2 , n = 1,2,3... 

(3.16) 

The potential V^(ar) is the rationally extended trigonometric Scarf potential, studied earlier in 
ref. [3]. The bound state solutions can be written in terms of exceptional X\ Jacobi polynomials 

[1. 2] 

^ n (x)-M n 2a-l-2/3smkx ^ {SmkXh n " 1 ' 2 '- 

(3.17) 

(c) m = 2 : In this case we have the following potential, energies and bound state wave functions 
involving X 2 exceptional Jacobi orthogonal polynomials 

V^ 2 \x) = k 2 [a(a — 1) + (3 2 } sec 2 kx — k 2 j3(2a — 1) sec kx tan kx 

4fc 2 [3(2a - 1)(2(3 + 1) sin kx - 2(3(2(3 + 1) - 8a(a - 1)] 

+ 2((3 + 1)(2P + 1) sin 2 kx + 2(2(3 + l)(2a - 1) sinfcx + Aa(a - 1) - 2(3 - 1 

(3.18) 

8(2/3 + l) 2 k 2 cos 2 fcx[2(l + /?) sin kx - 2a + l] 2 



[2(/3 + l)(2/5 + 1) sin 2 kx + 2(2/3 + l)(2a - 1) sin kx + Aa(a - 1) - 2(3 - 1] : 

£^ 2) = ^(n + a-2) 2 , n = 2,3,4.. . 
and 

^)(x)=J^m (l-smkx)^\l + smkx)^ 



-8k 2 



2((3 + 1)(2(3 + 1) sin 2 kx + 2(2(3 + l)(2a - 1) sin kx + Aa(a - 1) - 2(3 - 1 

x p( a -M^M-2) (sinA;x); n = 2)3j __ > 



(3.19) 



respectively. 



Case II: c < 0, extended hyperbolic Scarf potential family 

Let c = —k 2 , k 7^ 0. In this case we have g(x) = ismhkx with d = 0. The corresponding 
potential (we rename them as U^ m \x)) reduces to 

Jj{m) ^ _ _ k (2a +2b -1) ^ _|_ ^M^Lzili ^ x ^ _|_ 2k 2 17l(a — b — 171 + 1) 

,/2^ /. > iNF L u , c 7 , i \ • ■ 17 l p m-i\i sinh fcg) 
+ /e (a — o — m + 1) [a + b + (a — b + l)i sinh /cxj — (_ a-1 



P m ' (i sinh kx) 



k 2 (a — b — m + l) 2 cosh 2 fcx 



Pm-i^ (i sinh fcx) 



2 



— OO < X < oo 



-Pm a 1,6 (« sinh fcx) 

(3.20) 

The bound state wave functions and energy spectrum of the Schrodinger equation for this potential 
U (jn \x) are given by 

e>w = ^' (i -' s ""f i !:/; + ;t'; fa)l+ ' ^e****), (3.21a) 

Pm ' (2 smh foe) 

-)■■ -^W,-^-^) (3,,, 

respectively. The potentials U^ m '(x) are real provided a = b. Like the trigonometric Scarf potential 
family, the potentials U^ m \x) given in (3.20) are infinite in number and can be interpreted as the 
rationally extended polynomial-dependent hyperbolic Scarf potential family. Since the potential 
is free from singularity in — oo < x < oo so the wave function should vanish asymptotically 
at ±00. This imposes an upper bound on the quantum number n < (m — a+ ^ +1 ) . Also all 
the rationally extended hyperbolic Scarf potentials are isospectral with energy spectrum E^J 1 ' = 
-f {2n' + a + b+ l) 2 , n' = 0,1,2... 



C. Infinite family of non-Hermitian Hamiltonians 



In the previous section III B, we have considered the constant of integration d to be zero. The 
nonzero value of d corresponds to the shift of the coordinate x. In the Hermitian case this co- 
ordinate shift is not so important as they do not influence the energy eigenfunctions and eigenvalues. 
However, in PT-symmetric quantum mechanics imaginary coordinate shift plays a significant role. 
This gives rise to complex potentials with entirely real spectrum. Here we shall examine how the 
potentials V^ m '{x), U^ m \x) and associated energy spectra behave if we allow purely imaginary 
value of d. For this, we set d = ie and c = k 2 , in which case the equations (3.9) and (3.10a) 
reduce to 

g(x) = sin(kx + ie) (3.22) 
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and 

7' m '(i) = k <y2a — sec 2 (A;x + ie) — k ^ 2 ~ a ^ sec(kx + ie) tan{kx + ie) — 2k 2 m(a — b — m + 1) 
— k 2 (a — b — m + l)[a + b + (a — 6+1) sin(kx + ie)} 



PL a 1,6 1 ' ) (sin(/ca; + ie) 



^ k 2 (a — b — m + l) 2 cos 2 {kx + ie) 



P^_°i^(sm(A;x + ie)) 



oo < x < oo 



_Pir a_1 ' 6_1) (sin(A;x + ze)) 
£^ m) = ^(2n - 2m + a + b + l) 2 , n = m, m + 1, m + 2, ... 

(3.23) 

Like the Hermitian case, the potentials V^ m \x), m = 0,1,2,.. are infinite in number and all the 
members are isospectral. Here, a, b are chosen to be real otherwise the solution of this potential 
will contain Jacobi polynomials with complex indices and complex arguments. Such complex 
polynomials are not suitable for physical applications because of their non-trivial orthogonality 
properties which depend on the interplay between integration contour and parameter values. It 
is to be noted here that the non-Hermitian Hamiltonians with the potentials V^ m \x) have 
entirely real bound state energies This can be proved in the following way. Let us define a 

Hermitian, positive definite operator p as 

p = e ip p = _ijL (3.24) 
ax 

which has the following properties [44] 

pxp _1 = x-p pvp - l = V) pf(x)p~ l = f (^x - j). (3.25) 

With the help of the operator p, we have the following similarity transformation 

pV {m) (x)p~ 1 = V {m) (x) (3.26) 

This shows that the non-Hermitian Hamiltonians are quasi-Hermitian with respect to the 

Hermitian, positive definite operator p. The equivalent Hermitian counterparts are the Hamiltonians 
with the potentials V^ m '(x) (3.10a) corresponding to d = e = 0. The energy eigenfunctions of 
the potentials V { - m \x) can be obtained using ip^ = p~ lr ^n{%) = 4>n {x + ie), where i/in(x) 
are given by equation (3.11). It will not be difficult to show that the Hamiltonian is also 

pseudo-Hermitian i.e. 

rjH^r]- 1 = (3.27) 
with respect to the positive definite operator 

rj = p 2 = e f p . (3.28) 
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PT-symmetry of the potential V^ m '{x) cannot be achieved, in general, if d has a non-zero real 
component, because the finite shift along the coordinate x makes the potential different from 
its VT counterpart. Now, we shall try to determine under what condition the above potential 
remains VT invariant. A necessary requirement for a Hamiltonian to be PT-symmetric is that the 
potential should satisfy following condition 

V^ m) *{-x) =V {m \x). (3.29) 

This implies that the potential V^ m '(x) is PT-symmetric iff a = ±b. Hence, we have derived 
infinitely many non-Hermitian rationally extended trigonometric Scraf potentials with entirely real 
spectra and whose solutions are associated with exceptional X m Jacobi EOP. The non Hermitian 
variant of the hyperbolic Scarf potential family can be obtained in a similar manner and is given 
by 



!/' m '(i) = — k ^ 2a +2b — — sech 2 (/;;x + ie) + i k ^ b 2 a - sech(fcx + ie) tanh(/cx + ie) + 2k 2 m[a — b 
+ k 2 (a — b — m + l)[a + b + (a — b + l)i sinh(A;x + ie)} 



Pm-x* sinh(/cx + ie)) 
Pm. a 1,6 (i sinh(/cx + ie)) 



k 2 (a — b — m + l) 2 cosh 2 (kx + ie) I {i smh{kx + ie)) 



2 lpi- a - 1 ' b - 1) (^sinh(fcx + ^e)) 



— OO < X < oo 



(3.30) 

These complex potentials U^{x) are also quasi- and pseudo-Hermitian with respect to the oper- 
ators p and 7] respectively which are defined in eqns.(3.24) and (3.28). Therefore the potentials 
[A m )(x) have entirely real energy spectrum as of its equivalent Hermitian analog U^ m '(x) given in 
equation (3.20). Moreover, the complex potentials U (jn \x) satisfy the PT-symmetric condition 
U( m >(-x) = U^(x) for all real values of a and b. 

IV. SUPERSYMMETRIC SHAPE INVARIANCE APPROACH 

Here, we combine the results, obtained in section III, with supersymmetric quantum mechanics. 
Supersymmetry makes use of two linear first-order differential operators A^ mS)± = q=j- + W^ m \ 



where W^ m '(x) is the superpotential and generally defined [50] in terms of the ground state wave 
function ipm (x) as 

^\x)^exp(- J W [m \r)d^j (4.1) 
The two operators A^'^ give rise to two partner Hamiltonians iy( m )± 

H (m)T = A (m)± A (m) T _ ^_ + y{m) T _ (43) 

ax 2 
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where e is the factorization energy and V^™^ are the two partner potentials given in terms of 
superpotential 

V {m)± (x) = W {m \x) 2 ± W {m) \x). (4.3) 



The partner Hamiltonians have the same spectrum except the zero energy state i.e. E ( ^ l 



E 



and E, 



(m)- 



0. If the zero energy state of one of the Hamiltonians H^™^ is known, 
the eigenfunctions of the other Hamiltonian can be determined using ~ A~ip^{~ and 

V ; n ?1 ' ) ~ ~ A + ip^^~. The two partner potentials V^ ± (x) are said to be shape invariant if they 
satisfy [51] 

V M+ ( X; ao ) = V M- ^ + R ^ ( 4 4) 

where a x = f(a ) and R(a ) is independent of x. The beauty of shape invariance property is 
that whenever two supersymmetric partner potentials are related by shape invariance condition, 
the energy eigenvalues and the eigenvectors can be determined algebraically [50]. 



For convenience we identify the potential V^ m '{x), given in equation (3.10a), with 
The ground state solution [follow the equation (3.11)] of this potential can be written as 



m)- 



[X 



4 



m)- 



(1 — sin kx) 2 + 4 (1 + sin kx) 2+4 



1 - 



P, 



(-0-1,6) 



m—l 



sin kx) 



p(-a-l,6-l) , 



(4.5) 



(sin kx) 

Using eqns. (4.1), (4.5) and after some algebraic manipulations with the help of the recurrence 
relations (2.6) we obtain the superpotential as 



W^(x) = k -^seckx+ k{a+ ^ +1) 
k(a — b — m + 1) cos kx 



tan kx 

P^-i ] '(sin kx) 



P, 



(-0-1,6+1) 



m— 1 



sin kx) 



P 



(-0-1,6-1) , 



sin /cx) 



,(-0-2,6) , 



sin kx) 



(4.6) 



From equations (4.3) and (4.6) we get the simplified expressions of the partner potentials V^ T 



as 



y( m ) (a, b, x) = k ( - 2a +2b — — sec 2 kx — k ^ 2 a - sec kx tan kx — 2k 2 m(a — b — m + 1) 

-Pi^i' 6 ' 1 (sin /cx) 



— A; (a — b — m + 1) [a + 6 + (a — b + 1) sin kx] 



P 



(-a-1,6-1) 



fc 2 (a — b — m + l) 2 cos 2 kx ( P^l b \smkx) 



sin kx) 



(4.7) 



fc 2 (a + 6+l)' 



p(-a-l,6-l) 



sin fcx) 



and 

V (m)+ {a,b,x) 



k 2 [2(a + l) 2 + 2(6+l) 2 -l] 



sec kx 



k 2 [(6+l) 2 -(a+l) 2 ] 



- k 2 (a - b - m + 1) [a + b + 2 + (a - 6 + 1) sin kx] 



P 



-a-l,6+?)/ . 



sec kx tan /cx 



in— 1 



sin fcx) 



(-a-2,b) 



+ 



fc 2 (a — 6 — m + l) 2 cos 2 kx ( P^-i 



(-0-1,6+1) 



sin kx) 



sin fcx) 



,(-0-2,6), 



sin kx) 



2fc 2 m(a - 6 - m + 1) * ^ 



(4.8) 
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The potential y( m )~ matches, apart from an additive factorization energy e = k ( a + b+1 ) i with the 
potential (3.10a) obtained in section III. The potential V^ m \x) is another infinite set of exactly 
solvable potential. It is not very difficult to show that the bound state solutions of this potentials 
are also associated with X m Jacobi Polynomials. Also, it is very easy to check that the two partner 
potentials V^ T given in equations (4.7) and (4.8) are connected to each other by 

V (m)+ (a,b,x) = V (m) -(a + l,b+l,x) + k 2 (a + b + 2) (4.9) 

i.e. they have the translational shape invariance symmetry. Analogously one can study other 
generalized potential families, obtained in section III, in the framework of supersymmetric quantum 
mechanics to show that they have shape invariance symmetry. 

V. SUMMARY AND OUTLOOK 

We have obtained infinitely many exactly solvable Hermitian as well as non-Hermitian trigonomet- 
ric Scarf potentials V (ja \x) and hyperbolic Scarf potentials U^ m '(x), m = 0,1,2.... The bound 
state wave functions of these potentials are associated with the exceptional X m Jacobi polyno- 
mials. The potentials belonging to a particular family V^ m '(x), U^{x),V^ m '{x) or U^ix) are 
isospectral to each other. The supersymmetric partners of the potentials belonging to the above 
families possess shape invariant property. The non-Hermitian Hamiltonians involving the complex 
potentials V^ m '{x) and U ijn \x) are shown to be quasi-Hermitian with respect to a positive definite 
operator p = e^ p which affects the imaginary shift of co-ordinate x — > x — lp. This implies that 
corresponding energy spectra are all real. It has also been shown that for real a = b the potentials 
V^ m \x) are PT-sym metric whereas the potentials U^ m \x) are "PT- symmetric for all a, b. 

It is well known that the non-Hermitian conventional hyperbolic Scarf potential has both 
broken as well as unbroken PT-symmetry depending on the values of potential parameters [40]. 
Therefore it will be interesting to study whether there exist same type of PT phase transition in the 
generalized non-Hermitian hyperbolic Scarf potentials obtained here. Also the scattering state of 
these potentials are worth investigating to search for the existence of spectral singularity [52, 53]. 
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